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ABSTRACT 


A cylindrical acoustic resonator which has been developed by Dr. 
William J. Toulis of the Signal Propagation Division of the United States 
Navy Electronic Laboratory for measurement of the velocity and absorption 
characteristics of an enclosed fluid medium is discussed. The cavity is 
constructed of compliant tubes, thus presenting a pressure release surface 
at the boundaries. Theory is developed indicating that the resonant fre- 
quency of the cavity is directly proportional to the velocity of sound 
in the medium while the Q of the resonator can be used to determine 
attenuation. 

The writer wishes to express his appreciation for the guidance and 
encouragement of Professor Herman Medwin of the United States Naval Post- 
graduate School. He further wishes to acknowledge the assistance given 
him by Dr. Toulis and other members of the Signal Propagation Division 
of the United States Navy Electronics Laboratory in this investigation. 
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TABLE OF SYMBOLS 


a - cylinder radius 


A - vector potential 


B - bulk modulus of elasticity 


C,Co - speed of sound, reference speed 


Cp,Cy - specific heats at constant pressure, constant volume 


dj - rate of deformation tensor 


D - material differential operator 


En, Ey,E, - total energy, internal energy, kinetic energy 


frequency 


extraneous force 
NA vo 


-P C4? a frequency dependent ratio 
o Lo 


J NCp 
designation of imaginary axis 


» none-frequency dependent ratio 


phase constant - _ 
cylinder length 
number of pressure antinodes in the axial direction 


viscosity number = \\ = 2 + hs 


P; Pl»P2 ~- incremental, or sound, pressures; first and second order 


q 
a 
Q 
Ss 


pressures 
heat flux vector 
Stokes radiation coeffictent (see pp. 8) 
quality factor 
clesed surface 
specific entropy 
absolute temperature 


time 





ti, ° stress tensor 
Uiy - residual stress 
V - volume 
V - particle velocity 
"— - specific volume 
W - energy stored 
of - attenuation constant 
B - coefficient of thermal expansion 
‘E - ratio of specific heats - Cp/Cy 
6a4- Kronecker delta 
& - specific internal energy per unit mass 
6, - first-order variational component of temperature 
KK - thermal conductivity 
\ - dilatational viscosity coefficient 
i - dynamic shear viscosity coefficient 
~P 3 Pip2- total density; first and second order variational components 
CJ - time dependence variable (O #xAX + W) 
D - scalar velocity potential 
Dy- viscous dissipation function 
w - angular frequency 
Say - nth zero of Bessel function of order m 
V ©) - gradient 
V -© @- divergence 


V7 XO - curl 


vi 








1. Introduction. 

The success of sonic echoeranging methods depends, in part, upon the 
velocity and attenuation profiles in the transmitting medium between the 
detecting ship and the target. The evelopment of an instrument which 
could directly measure these properties of the sea would be of great 
assistance. 

The velocity of sound in sea water under varying physical conditions 
has been investigated both experimentally and theoretically by many pere 
sons. For the experimental determinations, much of the early work was 
pointed toward accurate determination of depth by sonic sounding methods. 
Transmission time between a transmitting and receiving hydrophone was 
measured in an area where salinity and temperature were as uniform as 
possible, and in deep enough water so that bottom reflections were mini- 
mized, In this type of determination, measurement of the distance is 
the limiting factor. 

More recently, a number of velocity meters has been developed for 
this purpose. The most successful for laboratery use has been the inter- 
ferometer, where a standing wave is set up in a cell by utilizing a 
crystal and reflector or two reflectors. This then provides a measure of 
the wave length and hence, the velocity. However, no completely satis- 
factory operational-type system is presently available.! The most widely 
used device to determine sound velocity profiles is the bathythermograph, 
which measures temperatures rather than velocities. Its use depends on 
the fact that the velocity of sound is largely determined by the tempera- 
ture of the sea, and only to a lesser extent by salinity and pressure. 
Urick? obtained good general agreement between velocity-depth plots 
obtained with an interferometer and those computed from bathythermograph 


and thermocouple observations for sea states of zero to three and depths 


Lu. Sussman, Study and Evaluation of ONR-Sponsored Sound-Velocity 
Meters, USL Research Report No. 299, 1 March 1956 


2R. J. Urick, An Acoustic Interferometer For the Measurement of 
Sound Velocity In the Ocean, USNR & S Lab. Report No. 5-18 


as great as 150 feet. However, the bathythermograph is insensitive to 
temperature microstructures, 2 and as such, does not give an accurate 
indication of transmission characteristics. 

A great deal of continuing effort has been given to the theoretical 
determination of the velocity of sound from the physical constants of 
the medium as given by the theory of thermodynamics .°»% 7 Much of the 
work is, of necessity, empirical in nature and has been complicated by a 
lack of sufficiently accurate experimental data.8 However, recent im- 
provements in laboratory measuring techniques should rapidly improve this 
situation. ? 

Techniques for measuring attenuation of sound by sea water are even 
less satisfactory. The absorption at low frequencies is sufficiently 
small to make measurement difficult without considering a large expanse 
of sea. Here again, determination of distance is a limiting factor. 
Further isolation of the specific causes of attenuation is next to ime 
possible without a device which can determine the absorption character- 
istics of small samples. 

The velocity meter under discussion depends, for its operation, on 


the fact that the resonant frequency of a water column in a particular 


3Urick, ibid. 


“5 Liebermann, The Effect of Temperature Inhomogeneities in the 
Ocean On the Propagation of Sound, J. Ac. Soc. Am., 23, pp 563 (1951) 


5s, Kuwahara, The Velocity of Sound in Sea Water and Calculation 
of the Velocity for Use in Sonic Sounding, Japanese J. Astron. & Geophys., 
Vol. 16, No. 1 (1938) 


Op, J. Matthews, Tables of the Velocity of Sound in Pure Water and 
Sea Water for Use in Echo-Sounding & Sound-Ranging, Hydrographic Dept. 
Admiralty, London, H.D. 282 (1927) 


V. A. Del Grosso, The Velocity of Sound in Sea Water at Zero 
Depth, NRL Report 4002, 11 June 1952 


SR. T. Beyer, Formulas for Sound Velocity in Sea Water, J. Mar. 
Research, Vol. 13, pp 113 (1954) 


SV. A. Del Grosso, E. J. Smura and P. T. Fougere, Accuracy of 
Ultrasonic Interferometer Velocity Determinations, NRL Report 4439, 
6 December 1954 





mode of vibration is directly proportional to the velocity of sound in 

the water column. The technique employed is not new in itself. A hard 

wall cavity resonator was developed by the University of Michigan in 1952.19 
However, this device was extremely sensitive to wall cleanliness due to 

the relatively high sonic pressures at the boundaries, characteristic of 
rigid walled devices. Thus, although the instrument provides an accurate 
measure of the velocity when used in the laboratory, it is very difficult 

to use under operational conditions. 

Use of a pressure release walled cylinder suggests itself; and, in 
fact, thin walled cylinders have been used for measurement of the atten- 
uation of marine sediments.41+12 Recent developments in the field of 
compliant structures by Dr. William J. Toulis of the Signal Propagation 
Division of the United States Navy Electronics Laboratory, San Diego, 
California, enabled him to adapt this technique to the present operational 
device, where, by measuring the resonant frequency and the '"Q" of the 
cavity, the desired characteristics of a small sample of water can be 
determined to a high degree of accuracy. 

The purpose of this paper is to discuss the theoretical aspects of 
the pressure release walled cylindrical cavity constructed by Toulis, 


including a somewhat systematic study of the loss problem. 


10p. kK. Brown, The Development of an Underwater Sound Velocity Meter 
Using a Cylindrical Resonator, Proj. 2021, final report, 30 April 1953, 
Eng. Research Inst. Univ. of Michigan 


llw. J. Toulis, Theory of a Resonance Method to Measure the Acoustic 
Properties of Sediments, NEL Report 676, 5 March 1956 


I2, Shumway, A Resonant Chamber Method for Sound Velocity and 
Attenuation Measurement in Sediments, Geophys. Vol XXI, No. 2, pp 305- 
319, April 1956 


3h 





2. General Wave Equation. 

Before proceeding to the actual analysis of the velocity meter, it 
is necessary to develop the equations governing the transmission of 
acoustic waves through a fluid medium. It is assumed that the medium is 
locally homogeneous, continuous except at the boundaries, and isotropic 
when at rest. It is further assumed that the condition of the medium can 
be fully described by establishing the functional dependence of pressure, 
vector velocity, density and temperature on time and the space coordin- 
ates. A set of equations consisting of 

ea) the continuity equation expressing the conservation of mass 

b) a force equation based on conservation of momentum 

c) a heat exchange equation based on conservation of energy 
will be developed in their cemplete form to present a rigorous point of 
departure for making the necessary approximations for the study of small- 
signal effects. 

To derive the first of the required relationships, consider an 
elemental volume V bounded by a closed surface S$. The amount of fluid 
crossing the surface in time AU is given by the surface integral: 

outflow « av), PVedS 


where 


V «= particle velocity 
- = point fluid density 


During this same period, the quantity of fluid originally contained in V 
will have diminished by the amount: 


loss = ats 3 dV 


where the negative sign is selected because density is a decreasing 
function of time. Since mass must be conserved, net outflow must equal 
the loss of mass and the two expressions can be equated. The surface 


integral is then converted by use of the divergence theorem giving: 


J [v-ev) + 2£]= 0 


Since the integrand is continuous and the volume arbitrary, this may be 


written as: 


>? : oe 

an ye {(-—P\/) = 

oe (we O 2.1 
which is one form of the well-known equation of continuity. 
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Considering the density is a function of both space and time, the 
total derivitive of f with respect to ( can be written in vector 


notation as: 


De Vaan es 
Bee Te + VV P 2 


13 i> 


utilizing the notation of Stokes”~ where Dt denotes the material de- 


rivitive and is defined by 
De 2 ) ss We grad ( ). 


Dt a 6 
Substituting from equation 2.1 into 2.2 gives 
Df U7 
= 4+ PTV = O 2.3 
Dae 


which is Euler's form of the continuity equation. 
To develop the second of the basic relationships, consider Cauchy's 


first law of mot fens” ° 
Divan ye 4+ ODAd 
P(Bus | = PFs OK 4 ove 


where: 
FF, = extraneous force vector, a function of position and 
time | 
Va = velocity component of 
Uae stress tensor. 
For an isotropic mediwm where stress is a linear function of the rate of 
deformation, the stress tensor can be resolved into a sum of a scalar 


and a residual stress tensor as follows: 
Cay = — pBay t+ Uar CBisy = 16 7s) 2.5 
where 
JP = pressure 
LAAse residual stress 


Dag = Kronecker delta 
POA = hydrostatic tension. 


13¢, Truesdell, The Mechanical Foundations of Elasticity and Fluid 
Dynamics, J. Rat. Mech. & Anal., 1, pp 125, (1952) 


l4r. vy, Hunt, Notes on the Exact Equations Governing The Propaga- 
tion of Sound in Fluids, J. Acous. Soc. of Am., 27, pp 1022, Nov. 1955. 








To a first approximation, the residual stress can be expressed by 


the linear terms of a power expansion in the viscosity coefficients as: 


Urge Ndyp Gag + 2p das 2.6 
where: 


= rate of deformation tensor 


dynamic shear viscosity coefficient 


dag 
B 
d # dilatational viscosity coefficient 
and a dM is defined by 


OYA ov ) 5 
Jag = + (BY Y aa a 
Combining the equations 2.5, 2.6, and 2.7 and substituting into 2.4 yields 
the vector force equation: 


FQV)= fF - F(V-V)V - VP +(Ateyp)v(V-V) 
~UIX(VXV) + (7° V)VA + 2(VN-V)V 2.8 
+VNX(VXV)- 


To develop the third of the fundamental equations, it is necessary 





to establish a series of relationships based on the law of conservation 
of energy. The total energy, E+, of a given system can be expressed 4s 


the sum of the kinetic energy, Fk , and the internal energy, E+, where 


Er=\ fe dU a 
‘U - specific volume = '/-p 
€, - internal energy density 
\V - volume of the system. 
The material rate of change of the total energy is then the sum of 


the rates at which work is done upon the system; that is: 
O(Ex+Es) _ 
Dt =AEK SEs) =J Pe AVA + btrgvad Sa - pa 2.10 


where 
q = heat flux vector. 

This postulates that the rate of doing mechanical work can be divided 
into two portions: chat arising from the extraneous force and that aris- 
ing from the forces of material continuity expressed in terms of an 
equivalent stress vector acting on the surface 5S. The third term ex- 


presses the rate at which thermal energy enters or leaves the volume. 





Applying the divergence theorem and utilizing equation 2.4 reduces 2.10 
to 


Due = 
PURE) = Vag das OAs 2a 


This is the Fourfer-Kirchhoff-C. Neumann energy equation. Using the 


relationships from equations 2.5 and 2.7, 2.11 can be written as 


DE ODA “JV = 
== 4+ PAA + PHVOV=a U Dale 
se. See Se: 
Considering the internal energy density, © , to be a function of S and 
U gives 
ae = (2&) (23) re (2& \ (ov 2.13 
jhe aS /y\ Dt O@V/S \DUB 


From the first law of thermodynanics, 


26. a (25) z 
(26 S P aB/y a, 
Thus, equation 2.13 can be written as 
we. s DSL Wy 
c=. = - TB eV 2.14 
Defining ®,, , a viscous dissipation function as 


by, = Arr d Ax 


and eliminating PRs between equations 2.12 and 2.14 gives: ° 
Ds Oe 215 
aianie Dt $., A 


Next, the effects of radiation and conduction must be considered. 
Since the divergence of Q gives the energy transferred away from the 
volume element, it must account for energy loss by either conduction or 
radiation. 


For conduction, the Fourier relation is used: 
bas Oaaie 
(24) cond a hy SOR 


A- thermal conductivity. 
Taking the divergence of each side gives: 
a 
2(arcond_ — as T-(WT)-(VA) 2.16 
OKA 
For the radiation effects, Newton's law of cooling states: 


2(Ar)vad = pCya(T- T) 2.17 


oO XA 


where 





where 
- radiation coefficient introduced by Stokes, Phil. Mag. 


(4) 1, pp 305-317 (1851) 
Cy - specific heat at constant volume = TAs), 
(‘T - Te) - local temperature excess. 
Before proceeding, the term T(2S) appearing in equation 2.15 must 
be expressed in terms of VY ,V , and J’. Following the methods of 
Hunt !? , this is given by 


eT (82) = #Cv(r-) ZY + BT 2 
TS Ce/cy 
Cp=a ae 


BP - coefficient of termal expansion 


PCYAT) 5 
Combining equation 2.18 with 2.15, 2.16 and 2.17 gives the third of the 


required cas 


acy 





Vv 2 Ma a v" 
aa eed H] vk P 


-V TVA + POE (T- eae = O 


ISnunt, op. cit., pp 1023. 








3. Small Signal Acoustic Equations. 

The small signal acoustic equations are obtained by arbitrarily 
'linearizing'’ the general equations of motion. For this purpose, the 
dependent variables will be replaced by the sum of their equilibrium or 


zero-order values and the first-order variational components. Thus, 
P= fo +f P= Pot p 
= lope Vi = Orre 


On the assumption that ~~, , ars and Po are constants and that \/ 
is itself a first-order term (i.e. Y,= © ) then: 


Ge ver V*?P - VY 
Vee Sy Ven= Nig, 
Substituting these relations into equations 2.1, 2.8 and 2.19 and re- 


taining only the first-order terms gives: 


a te ae =O = 
(4) + Vi = (WN) V(7-™) + NV X(VXV, ) = 3.2 
0 Cy een 4 Ez Cv (v- )/p| (VV) “KV. 3.3 


os mye O, = 0, 
where the extraneous force has been omitted since it primarily affects 


the equilibrium configuration and 


N= 2 + Ay » the viscosity number. 
Considering “> as a function of Be and de gives: 





| = oP > oF 3.4 
Since: 2 
te ee. 
— | 
and 


2+ Peau yay) 
ca, “_ he: ah = PP, 


equation 3.4 becomes 


ae 
VP, = V+ Tt fs V2 








and 3.2 can be rewritten as 
—— 2. 
2@8) +(FAfire-niggi]ya Nv) ,. 
+ OS; | =e Ox 

These equations (3.1, 3.3, 3.5) take into account the first-order effects 
of shear and dilatational viscosity and of heat conduction and their 
solution is discussed in Section 5. They can be reduced to the classical 
small-signal wave equation by neglecting the effects of heat conduction 
end viscosity (K= 9 «AJ =) =Q) and assuming adiabatic behaviour of 
the medium giving 

op 

ae 





fs 
= Gay 4D 326 


10 





4. Solution of the Simple Wave Equation. 

Before considering the more complex equations derived in the pre- 
ceeding section, a solution will first be obtained for the small-signal 
wave equation. The results obtained will serve as a valuable guide when 
interpreting the solutions obtained in later sections. 


The simple wave equation is given by: 


oe a ee Al 
ude acme | 


Pe 
Where \f is the Laplacian operator and C is the ErOCMARAS® sound 
AT 
velocity. Assuming a time dependence of the form € where T=Q+AQM, 
equation 4.1 takes the form of a scalar Helmholtz equation: 


(V+ K*)-p alee ae 
where : 
2s. JOE ee a ON 
Nom CG ep eK NO 


Using the method of separation of variables with the Laplacian operator 


in cylindrical coordinates gives: 
i) 


R" Rg | & ' 4 oy _ 
Bm tetre te +h =° 4&8 
where rR" denotes as » denotes =, ’ 


for the method of separation of variables, equation 4.3 may be rewritten 
rf \ 1 
R ee | & Le Ea 
ae +R Tee +k = De 
This et produces two equations: a 
\ { 
na See eK 


2 
Zo + Kz Z= O 4.4b 


where equation 4.4b is the familiar differential equation of simple 


Qu 


etc. By the standard argument 





N 
~ 


4.4a 


harmonic motion with the solution: 


cS : A Ka @ \ 
In general, Ke can be considered complex; that is 
1 { 
Ka= AK + 4 


Thus, equation 4.5 may be rewritten as 


aie 4.5 


ee \ -%>, = 
x. = (Ax, 628 Ka Zz + PBK. 2M Ko. ze) = 45a 


11 





Equation 4.4a may be separated utilizing the same techniques giving: 


@' + m~@ = QO 4.6a 

ae 9 Ce 

rR 2 as R pees (0h a= K Kz = © 4.6b 
where solutions to 4.6 may be written in the form: 

@# = Cm COS(M G+ Y ) 4.7 


Since oy must be single valued in © ,\Y\ must be an integer. 
Equation 4.6b is recognized as a form of Bessel's equation providing 


a solution (with Yy, as an integer) of the form: 


X= Dyn, Ke USC a En kr Nua (Kye ©) soe 


where ue (Kr)is a Bessel function of the aes as and Nyy (KY) is a 
Bessel function of the second kindl®, and Ke Ke ee The boundary 
conditions which must be satisfied by equations 4.5a, 4.7 and 4.8 are: 
a. The pressure is finite throughout the interior of the cylinder, 
b. The walls are assumed to be pressure release surfaces requiring 
that Jo=© at the boundaries of the cylinder. 
Selecting a coordinate system with 2+ © at the bottom of the cylinder, 
and for a cylinder of radius a and length |_ , the following conditions 
apply: 
a. Coefficient of the cosine term in equation 4.5a must be zero 
Since A> is zero at Z2=0. 
b. Periodicity in Z is given by K=20 where Q. = 1,2,3-- since /p 
is zero at Z#=(_. 
¢. The coefficient of the second term of equation 4.8 must be zero 
since Bessel functions of the second kind cannot be present with 
YP finite at the origin.) 
d. Since the selection of the coordinate system is purely arbitrary 
the ©=0 axis can be positioned so that the phase angle in 
equation 4.7 is zero (Y=Q) with no loss of generality. 


1OucLachlan, Bessel Functions for Engineers - 2nd Ed., Oxford at 
the Clarendon Press (1955) 


L7McLachlan, op. cit. 


12 





Thus, the solution of equation 4.2 takes the form: ot 
A 


= —-Xe 
areeeeen yn, Ka ce “SIn( Kez) COS (MA) Sin(KeHE 49 
The final condition to be satisfied is that the pressure must be zero 


at = O- Thus 


Jam (Ky @) = O 


Defining a quantity S wy such that 


Sinn = Ke 4.10 


then 


Seals) O 4.11 


S= Sm a 
where {\ is an integer denoting the nth zero of a Bessel function of order 


YY) . Further, it can be shown!® that all roots of equation 4.11 are 
real; that is, Syuy must be real. Thus 


Ka Kran 4.12 
v On 
is real. Since 
cm 2 ek 
Ke = K — Kg 
then 2. 
a rn 
2 2 
Substituting the complex forms for Ko and K gives: 
a 2 \ > x \ ay, 2 
—%> + Ko + PONS = Soe, Kk oe ZAK K- a 4.13 
Equating the real and —* panes of equation 4.13 gives: 
cx 1 an 
—Ae + i = —,* = K' ses . 4.14a 
\ \ 
Mp Ka = a, K 414d 


Combining equations 4.14a and 4.14b and solving for K gives: 


K= as Vee + (2p%)*- wT (Seay? At «7 | 4.15 


If &> is assumed negligable, 





kK = re + (20) 4.16 
18 


A. Gary and G. B. Mathews, A Treatise on Bessel Functions, Mac- 
millan & Co. Ltd. St. Martin's Street, London, pp 79, 1931. 


13 








* 
F ad 


where the negative root is rejected since only positive frequency is being 
considered. As might be expected, equation 4.16 is the equation for the 
wave number in a lossless medium. Substituting for kK" and Ka in equa- 


tion 4.16 gives: 





WY Qinn, aes (Atr)" 5 (Sms) 
Se ~ VEL Sa 


Therefore, fQyny> the resonant frequency of the cavity is 


ie s\(hy 1 [an-.)\* 4.17 


ny ON 
It is desirable to excite a mode of vibration in which the medium 


has no tangential components of velocity. For a medium having even a 
slight viscosity, such motion at the walls can introduce a much greater 
damping of free vibrations than that due to absorption of energy by the 
medium itself. Thus, only axially symmetric modes will be considered 
(i.e. m= 0) for the remainder of the analysis. For the lowest axially 
symmetric mode, Sry = So} =x «(0.7655) and Q = 1. This corresponds to 
a pressure antinode at the midpoint of the axis of the resonator. 

If attenuation is present, its effect will be to lower the resonant 


frequency slightly as indicated by equation 4.15. 


14 





5. Sound Absorption and Dispersion By the Medium, 
Investigation of the general mechanism of absorption by the medium 
will be based on a solution of the first-order equations obtained in 


section 3. They were 
af 


BU to (V7 Gun 
Po cv (82) +|fFoCvr-1)/2 |] (7-V) -Ky*e@ 
a oe O (3.3) 5.2 
Fo (2%) + (ae (Bg. vP—-vNviyv-v) 
(30) sei - Se 2) (3.5) ee) 


+HIX(7XV) =O 
where the subscripts indicating first-order effects have been omitted. 
Differentiating equation 5.3 with respect to time and substituting 
from 5.1 gor 2 gives 
oe - te&e V(V-V V) 
~ PNY (7-2 Or one ts Na Yo 


Considering that, hecoraie to Helmholtz's ae any vector can be 





5.4 


uniquely separated into the sum of a gradient of a scalar and the curl of 


a divergenceless vector, let 
Vs Vo +U0XA Spee) 


where 
d - velocity potential 
‘A - vector potential 
VA = 0 
Substituting equation 5.5 into 5.4 yields the following scalar and vector 


wave equations: 


e/ 2 
a ~B& Vo + So B 4,32 NY Ge) =0 5.68 


eA Py Ky XA = 6 5. 6b 


It is interesting to note that equation 5.6b which governs first-order 


vorticity is not a true vector wave equation but a vector analogue of the 


L9Morse and Feshback, Methods of Theoretical Physics, Vol. 1, McGraw- 
Hill Book Co., pp. 53, 1953. 


ae 








typical homogeneous diffusion equation. Since no source terms are present 
in 5.6b, this implies that, if first-order vorticity exists, it must be 
generated at the boundarles. This lack of source terms is due to the 
technique utilized in generating the first-order equations. Based on the 
assumption that the equilibrium values of the state variables are con- 
stant and that the viscosity coefficients are functions of these state 
variables, then it is consistant to assume that VNo eVAo =©. There- 
fore, the source terms involving the gradients of the viscosity coeffi- 
cients appear as second-order rather than first-order sources. This section 
deals specifically with absorption by the medium; thus, solution of 
equation 5.6b and its implications will be reserved for section 6. 
Neglecting the effects of radiation and assuming a time dependence 


AW 
of the form & » equations 5.2 and 5.6a ai 


W>[ Fete FAWN | + wX~Po b-Aw SB 4,8= Or Sosa7 


VO PoCvr-N/pl-KV"e + AWPLCVS = O 5.8 


Solving 5.7 for @ and substituting into 5.8, and introducing the shorte- 


NAW 


G = ~4 xz, a frequency dependent term 
Po Co 


¢ 


hand notation: 


and 
H «= eae which is non-frequency sensitive 
A) NCp ’ 
nae 
Ro Garce 


fo() +aG (itv H)] 9% — ES Gren - -at) |v 6: OF 


Equation 5.9 is biquadratic in the complex propagation constant and it is 


where 





GH = 


gives 


convenient to rewrite it as 
(v 7+ K,?) (w+ Ge) io 5.10 


The standard quadratic formula can be used to determine the reciprocal 


squares of the propagation constants as 


aA ® 2) 

K> -Co | _ 4A8- GH \ 

= ae Sie TH) STAG “o HYI + 5.11 
ia PLES” ae ( + ) MI [i+AG(i- TrAG FA] ‘) 


Kk,” 
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For the frequency range under consideration with water as the medium, G 
and H are small numbers and equation 5.11 can be expanded by power series 


methods to a 


al ra z13 + H(v- \]- 4 = Gli + H(e-1)] § 5. 12a 


K=O oon aT a | ai @ ed | Cee e] 5.12b 


Neglecting second-order terms gives: 


a oe ~ AINA 5.13a 
Ky 3 a {Ne + (6 Eee - 


OO Ae = HON 
Le (a feeee. [red (e-HU- = aE ae 5.13b 


Referring again to equation 5.10, we see that the solution will take the 





form of two waves referred to as types I and II whose wave numbers are 
respectively K, and Ke . Examination of equations 5.13 indicates that 
the two solutions might also be classified as "cempressional" and ''thermal" 
since the type I waves represent the solution to the ordinary wave equae- 
tion with a slight damping coefficient and the type II waves are generally 
the result of thermal conduction. 


From equation 5.13a the absorption and dispersion measures for the 


see 23) AN : eet) a K | 


(E = 1 tEWON/ oy 
c= Colt + SWHN/gRc)* | 


The attenuation thus determined is the Kirchhoff approximation for low 


type I waves are: 


x = 





frequencies. 

The wave length of the type II wave ( Ne am(2%, “) is very much 
shorter than that of the type I wave. Further, the type II waves are 
very rapidly attenuated at the frequencies under consideration with a 
decrement of Se aS , and the velocity of propagation through the medium 


is given by: 


17 








6. Wall Losses. 

To study the effects of wall losses on the system, consider the force 
equation derived in section 2 (2.9). The viscosity coefficients AU and » 
are assumed constant, and effects of the extraneous force, F, are assumed 
negligible. Further, since only small amplitude vibrations are to be 
considered, all terms involving squares of the velocity components are 
neglected. Equation 2.8 then becomes 


PY) = -9p s(\iw)V(V-V)FANV'Y 6.1 


Since it was shown in section 5 that the thermal waves are rapidly atten- 
uated, the effect of heat conduction at the boundary will be neglected. 
Thus, before proceeding further in the analysis, a direct relationship be- 
tween pressure and velocity must be established@9, Flow of mass out of 
any volume element reduces the pressure in a compressible fluid. In an 
elastic fluid, if compression is small, the rate of change of pressure is 


proportional to the divergence of the velocity; that is 


Ec s 
ee NE =l e 6.2 
SD BWV 
B - bulk modulus of elasticity 
Taking the time derivitive of 6.1 and introducing the expression given in 
equation 6.2 for @& ives: 
qd > & 


(EL) = vievv)ewr Ave) +yv 3h ss 


—— 


Again defining VV as 
Ws vb ~~ VYA 6.4 
and substituting equation 6.4 into 6.3 gives the two equations: 
$2 = -Nvx7xA 6.5a 
Fd 


ey / \ rat 
ape = OY Y tt ay a 
Assuming a time dependence of the form ate where G=W+AQ , equations 


6.5a and 6.5b take the form: 


20Morse and Feshbach, Methods of Theoretical Physics, Vol. 1, McGraw- 
Hill Book Co., pp. 162, 1953. 
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Vvx7xA + h* A =O Pie. 


(77+ K>)\iGe= ee 6.6b 
where K2 = ECGs 
ae om 


K2 = ot Cae 
Br aAT(ANtA 
Since ay satisfies a scalar Helmholz equation, the form of solution will 


be the same as equation 4.9, that is i 
. - ens a 
= SINKz= Jrn(Ker) yy) (wm 2) im 6.7 


where 


K 2 = om Co 
Ke = a eG Ca a 
Jn (Ke A) = O 


Equation 6.6a takes the general. form of a vector "wave'’ equation. Accord- 
q q 





ing to Stratton?! the three independent vector solutions to 6.,6a are given 
by: 
ae =V'¥ 6.8 
M=zaVx ay 6.9 
NF ree WO 6.10 


where O\ is any constant vector of unit length and NU, is a scalar funce 


tion satisfying the equation: 


Vy + h*y=O 6.11 
From 6.8, 6.9 and 6.10, 
7>M=7N =O 
ei = ey 
Since, by definition, Y A= O, 6.9 and 6.10 must be the complete solution 
of 6.6a; that is, a 


A= M +N 


21ly. A. Stratton, Electromagnetic Theory, McGraw-Hill Book Co., 
pp 393, 1941. 
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Equation 6.11 is again a scalar Helmholz equation and thus YU is 
given by: 


W = eO> Dz COG i= (gv) 


Salil) Sin 


Considering , the zero mode (m = QO) 


YW = Sty, DB UNG:a0) 6.12 
Substituting equation 6.12 into 6.9 and 6.10 and selecting QO. = Og gives: 
M= 7x Ga VY 
ee COS ee 6.13 
=. aa Sin oz J (av) 
N= w(vx®) colt 
VY) Wor gb oon mae Oe J (av) +42 g ei (ae) Sy > bal 


Taking the curl of 6.13 and 6.14 gives 


VXA= 0 jog 0 ba J, (ax) | 28 
+ bh Ag E aaa a = : e 


611 
“ a2lg* Solar) ae ve | 


From the symmetry of the problen: 
Vy (vr, 2) NO Co es %) 6.16a 


i 


i 


L 


Wes (7, #) = ~Va (Vv, L -2) 6.16b 
Thus, the rotational component of the vector velocity has the form 


VXA= Arlen g(bg sinba J, (gx) | 


at Ad [Deg (ne Cos Jy (gv) | 6.17 
+ AalCog (@? Tolar) os ba | 

b> 4 a” es \y * 

b> - “ae A= eae 


From an examination of equations 6.7 and 6.17, it can be seen that since 


K.= b , then AF K,. Thus Ju(Q cA) Soe 
The irrotational component of velocity is determined by taking the 
gradient of 6.7. Thus 


Ties Cele Ke (H Ke Si} Ke te ey (ue) | 6.18 
+ Or [Eg ve (Ke (05%, 2 Jo er) 
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A rational accounting of the actual wall losses is well beyond the 
scope of this discussion. Derivation of even a reasonable approximation 
for the pressure field existing at the boundary would be difficult and 
tedious. Even if these equations were known, the mechanism of acoustical 
energy transfer to, and the absorption by the compliant structure is not 
sufficiently understood to make application feasible. Thus, it is nec- 
essary to assume a simple model, considering that the solutions obtained 
will act as a guide in interpreting the actual behavior of the resonator. 

As a first approximation, the pressure release boundary condition 
is retained with the additional assumption that the cylinder is bounded 
by a smooth wall. It is further assumed that, for a viscous medium, there 
is no slipping between the medium and the wall. This implies the condi- 
tion that the tangential components of velocity will be zero at the 
boundary. 

Examination of equation 6.18 indicates that the irrotational com- 
ponents already satisfy this condition, i.e., 

Yr 7 O av Z-0,287b 
You= 6 at T= OO: 
Considering the rotational components, since Jiao)z 0, Vox must be zero 


so that 
Vqg2 0 AV Y= 


Further, Coq must equal zero for the longitudinal component, Vz , to 
satisfy the boundary condition. Thus, for this simple model, there is no 
rotational velocity component, and no firsteorder vorticity exists. There- 
fore, the solution obtained in section 4 is unaltered. 

A closer approximation to the actual conditions can be obtained by 
assuming that pressure approaches some small finite value rather than zero 
at the boundary (1.e., Jo( Keay #0 »>ankKbL =O). In this case, the 
coefficient Dog in equation 6.17 must again be equal to zero. However, 
from the two conditions: 


NEO Cue = OR 
Si N= CC ne = 


two equations are obtained: 


Zi 








E ke sin kK, H% TK) | + © log Sate bz a) (ar) = QO 6.19 
at 2=O and 2=L, and 


E|kzcos kz To(Keo) | es c| os ba Jo( ga) |= QO. 6.20 


Since equation 6.19 must hold for all values of \~ , it follows that 


Ka= = i, Equation 6.20 can then be rewritten as 


Pl 
Te (Ky = aa Jo (aa) 
= aa ae (Qo). 6.21 
Equation 6.21 again indicates that for perfect pressure release walls (1.e., 
ha) =O ) there is no firsteorder vorticity and no viscous loss of 
energy at the boundary. 

As the boundary deviates from the ideal case, the coefficient Cc’ in 
equation 6.21 increases accordingly and vorticity is generated. Thus, 
friction arises between the layer of fluid adhering to the walls and the 
immediately adjoining liquid in motion causing a dissipation of energy at 
the walls. 

This effect can be evaluated by considering that 

KA = Smun—- Emn 
where em is a small complex quantity which is dependent on e » Ree 
| ferring to equation 4.13, it can be seen that the addition of Caan will 
result in a decrease in the resonant frequency of the cavity as well as 
an increase in the total energy dissipation. 

A further implication of this analysis is that, since P #O atY=Q, 
a transmitted as well as a reflected wave exists at the boundary and some 


energy will be lost to the surrounding medium due to radiation. 
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7. Measurement of Attenuation, 

Before proceeding to apply the results obtained in the preceding sec- 
tions to the actual design of a cavity, some method of measuring the 
energy absorption by the enclosed medium must be determined. The quality 
factor or Q of a resonant cavity may be defined as: 

a 2x {energy stored in cavity) 7.4 
energy lost per cycle 
The mechanical energy stored in an incremental volume can be considered 


as the sum of the kinetic and potential energies; that is, 

E+= Ex + Ep dec 
where a fv" \/ 

Ep = Pi*/2.2c* dV 


Therefore, W, the total energy stored in the cavity, is determined by tak- 
ing the volume integral of equation 7.2; i.e., 
w= f Hedy ee an AV res 
Vv & Y Sogo 


For the present, wall losses are neglected and all attenuation is assumed 





ov 
to occur in the medium. Selecting a time dependence of the form e* where 
GQ is again equal to AX +O and substituting into equation 7.3 indi- 


—-2X%D 


cates that W decays at a rate determined by & where & is the atten- 


uation constant in nepers per second. Thus, the energy lost during one 


2% WW 


energy lost per cycle = ee 7.4 


period is given by: 


Substituting equations 7.3 and 7.4 into 7.1 gives 
2T7TfW ~~ _& 





2. Ok VW 2. A 
or 
Oo 
= nepers/second 15 


To convert & from a time to a space dependence, equation 7.5 is multi- 
plied by Ve: 


\ 





, co K . 
of 2cQ = ZQ nepers/unit length 7.6 


where K' is defined by equation 4.15 as 
ee 2 eae eee a ey, 2. ZS 
i 2 + (Bee )h— oF ( (Sean) 3 4 0%] 
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Substituting in equation 7.6 and solving for C% gives: 
he ee > 
A= Ka > 4+ (2xnn\* 





2Q 
29 \(t)* + @an* 1 


Thus, with the assumption that all energy dissipation occurs in the medium, 





the quality factor of the cavity can be utilized as a measure of the ate 
tenuation constant of the medium. 
In the case of an actual resonator, the effect of wall loss cannot 
be neglected, and the cavity is characterized by a total Q, Qy, given by 
pe rs 
CE een ee 


where G),, is due to dissipation by the medium and Qw is a correction 
due to wall loss. Qian » then, is the value which must be utilized in 
equation 7.7 to measure absorption by the mediun. 

Since the device under consideration will be operating at relatively 
low frequencies, determination of the correction factor, Qw > is greatly 
simplified. 

It is well knownthat, at frequencies below 10 k.c./sec, attenuation 
of sound waves in fresh water is considerably less than attenuation of 
sound in sea water. In fact, by comparison, attenuation by the fresh 
water can be neglected. Thus, if the cavity is immersed in distilled 
water and the total quality factor is determined, then 


—— 


a ~ ee 
Qt, Qwe 


to a close approximation. The subscript d is used to indicate measure- 
ment in distilled water. 


Q may also be defined as 
Q= i 
AF 





where 
{ + resonant frequency 
Afaf,-f, 


£,4,- frequencies of the upper and lower half power points. 


Therefore, 
eee (St \ _ = ta 
foe sae ci fg IW 


Considering that the resonant frequencies in fresh and salt water will be 


only slightly different, 
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a Ia ee 
Qyy - 
where 


At . band width between half power points for actual cavity 


in sea water 
AS - band width between half-power points in distilled water 
£ - resonant frequency in sea water. 


Substituting this value into equation 7.7 gives 


TAs = Se) 
: Boe 


and equation 7.8 can be used to experimentally determine attenuation by 
the medium. 





7.8 
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8. Cavity Design. 

In the design and construction of a soundevelocity and attenuation 
meter based on the principles discussed in the foregoing analysis, two 
general problems must be considered: 

a. Design of a cylindrical cavity presenting essentially a pressure- 

release boundary over the frequency range under consideration. 

b. Development of the "intelligence" system including a transducer 
to excite the water column to resonance, a receiving transducer 
to sample the generated pressure field and the associated 
electronic equipment for interpretation and display of the de- 
sired information. 

As a first step in the design, the general usage of the data obtained 
from the instrument must be considered to determine the desired resonant 
frequency of the cavity. Equation 4.17 can then be utilized to calculate 
the approximate dimensions of the cylinder. Care must be taken to select 
a length to diameter ratio that will minimize danger of mode skipping. 

Next, the general range ef sound velocities to be measured must be 
specified. This, in turn, determines how broad the Q of the cavity must 
be to present adequate sensitivity at the extremes of the operating range. 

With these parameters in hand, the formulas developed by Toulis22 
can be utilized to determine the nature and spacing of the compliant 
elements for construction of the cavity walls. This, then, outlines the 
steps to be carried out in the solution of the first portion of the 
design problen. 

In selecting and mounting the transducers for exciting the water 
column into resonance and sampling the generated pressure field, several 
factors must be considered. First of all, since the material generally 
used for enclosing the transducers is highly absorptive, it is desirable 
to mount them in a relatively low intensity sound field. Further, it 
is necessary to keep separation sufficient to minimize direct coupling. 

The design of a system for interpreting and presenting the data ob- 
tained is primarily dependent upon the actual operating requirements and 
will not be discussed. 


22, J. Toulis, Acoustic Refraction and Scattering With Compliant 
Elements, to be published. 
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9. Conclusions. 

Examination of the foregoing analysis indicates that the solution 
to the simple wave equation obtained in Section 4 describes, to a very 
close approximation, the conditions existing in an actual resonator. 
From Section 5, it can be seen that the effects of heat conduction on 
the resonant frequency and attenuation constant are negligible. Section 
7 indicates that, with walls closely approximating an ideal pressure- 
release surface, the effects of first-order vorticity can also be neg- 
lected. Thus, for design of a system employing the techniques under 
discussion, the simplified analysis presented in Section 4 should prove 
adequate. 

One particularly desirable feature of this system is that it mea- 
sures velocity and attenuation simultaneously. This should be of great 
assistance in determining the effect of air bubbles and other natural 
phenomena on sound velocity. 

Another desirable feature of this equipment is the frequency range 
employed. Since some doubt exists as to the validity of adapting sound 
velocity data measured at megacycle frequencies to much lower-frequency 
applications, an accurate device operating at these lower frequencies 


should be a welcome addition. 
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APPENDIX I 


CAVITY CONSTRUCTION AND TEST 

A cylindrical acoustic resonator sound-velocity meter was constructed 
by Toulis and preliminary tests were conducted by Toulis and the author 
from 6 January 1957 to 15 March 1957. This appendix outlines the methods 
of construction and instrumentation utilized and describes the test 
results. 

The pressure-release walls of the cavity under test were constructed 
by deforming commercially available aluminum tubing (Alcoa 3003-H-14, 1/2 
hard, 1'' 0.D., 0.036" wall thickness) into nearly flat or elliptic cross- 
section and welding the ends to prevent flooding. These tubes were then 
secured to a cylindrical frame with 2-1/2 inches spacing between tubes 
forming a closed cylinder. 

The water column was excited into resonance by a cylindrical barium 
titanate transducer mounted just below the compliant tube top of the 
cavity near the cylinder wall, and a similar transducer utilized to detect 
the pressure field was mownted directly opposite. 

For making velocity measurements, the output of the receiving trans- 
ducer was amplified and applied to the exciting transducer. A self- 
excited oscillating loop was thus fermed. The loop voltage was applied, 
through a frequency multiplier (X10), to a Potter 150 KC. Frequency 
Counter. 

For attenuation measurements, the amplitude of oscillation ts an 
indirect measure of Q. Direct determination was made by driving the cavity 
with a General Radio Beat Frequency Oscillator (Type 1304-A) and measure 
ing the bandwidth between the half-power points. 

During the tests, the output of the receiving transducer was applied 
to the amplifier through a calibrated potentiometer. The amplifier out- 
put was monitored and the gain (potentiometer) setting required to main- 
tain constant output was recorded as an indication of the attenuation. A 
block diagram of the system is shown in Figure l. 

To evaluate the ability of the instrument to measure actual sound 
velocity at sea, a number of tests were conducted in the coastal waters 
of San Diego, California. For comparison purposes and, particularly, to 


isolate the effects of aalinity, an additional set of runs was made at 
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the NEL testing facility at Sweetwater Lake. At each coastal station, 
the cavity was lowered by means of a bathythermograph winch, and the 
following data was recorded as the chamber was stopped at intervals on 
the ascent: 

a. Resonant frequency as indicated by the Potter counter. 

b. Temperature as determined by a thermistor bead attached to the 

side of the cavity. 

c. Gain setting. 

d. Depth as indicated by a wire metering device. 

The results of these tests’? indicate that the System tested pre- 
sents promise for development as a survey instrument for making sound- 


e 
fd 


velocity profiles. 


23 
J. E. Wille, Preliminary Test of A Cylindrical Acoustic Resonator 
Sound-Velocity Meter, NEL Tech. Memo., to be published. 
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